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Abstract 

We consider the lightest supersymmetric particle (LSP), neutralino in min- 
imal anomaly mediated supersymmetry breaking model (mAMSB) to be a 
possible candidate for weakly interacting massive particles (WIMP) or cold 
dark matter and investigate its direct and indirect detections. The supersym- 
metric parametric space for such a model is constrained by the WMAP results 
for relic densities. The spin independent and spin dependent scattering cross 
sections for dark matter off nucleon are thus constrained from the WMAP 
results. They are found to be within the allowed regions of different ongoing 
direct detection experiments. The annihilation of such dark matter candi- 
dates at the galactic centre produce different standard model particles such as 
gamma rays, neutrinos etc. In this work, we investigate the possible fluxes of 
such particles from galactic centre. The neutrino flux from the galactic centre 
and at different locations away from the galactic centre produced by WIMP 
annihilation in this model are also obtained for four types of dark matter halo 
profile. The possibility of detection of such neutrinos from galactic centre at 
the ANTARES under sea neutrino detector is also investigated. We have stud- 
ied signals from dark matter annihilations from different angles of observations 
for different spherically symmetric dark matter halo distribution models in the 
galaxy. We have compared our gamma ray flux results for four different halo 
models with the HESS experimental data. 



1 email: kamakshya. modak@saha.ac. in 
2 email: debasish . ma j umdar @ saha . ac . in 



1 



1 Introduction 



Cosmological observations like flattening of rotation curves of spiral galaxies [TJ, the 
gravitational microlensing [2J, observations on Virgo and Coma clusters [HIE], bullet 
clusters [5], etc. provide indications of existence of huge amount of non- luminous 
matter or dark matter (DM) in the universe. The Wilkinson Microwave Anisotropy 
Probe (WMAP) experiment j6] suggests that about 85% of the total matter content of 
the universe and almost 23% of the total content of the universe is dark matter. The 
general wisdom is that in order to account for the relic abundance of DM, a candidate 
for dark matter should be massive, very weakly interacting and non-relativistic (cold 
dark matter or CDM) particles. This allows the structure formation on large scales. 
In the present work, we consider such weakly interacting massive particles (WIMPs) 
[3 El El EED] to consist of the total DM content of the universe. 

The theory of Supersymmetry provides a solution to hierarchy problem and uni- 
fication of gauge coupling constants via renornalization group evolution (RGE). R- 
parity conserving SUSY also provides very naturally the lightest supersymmetric 
particles (LSP) to be a possible candidate for DM. In the present work, we consider 
neutralino to be the LSP and hence the candidate for dark matter. 

In SUSY models, R-parity, or some similar parity property, allows only an even 
number of supersymmetric partner particles to interact on a fundamental interaction 
vertex. This stabilizes the lightest supersymmetric particle (LSP), which becomes 
the cold dark matter candidate. Minimal Supersymmetric Standard model (MSSM) 
with softly broken supersymmetry is the main candidate of physics beyond the Stan- 
dard Model. Supersymmetry, if it exists, must be broken spontaneously. Dynamical 
or spontaneous breaking of supersymmetry at high scale leads to the soft Supersym- 
metry breaking terms appearing in MSSM as low energy remnants. This dynamical 
or spontaneous breaking is supposed to take place in some 'hidden' sector (HS) and 
this breaking is mediated to the observable sector (OS). This mediation mechanism 
leads to many interesting theories including gravity-mediation (SUGRA) with grav- 
itino mass (ms) (~ ITeV), Gauge mediation (GMSB) with mz < 1 TeV, anomaly 
mediation with ms ~ 100 TeV. The superconformal Anomaly Mediated Supersym- 
metry Breaking (AMSB) mechanism is one of the most well-known and attractive 
set-ups for supersymmetry breaking because, 

1. the soft supersymmetry (SUSY) breaking terms are completely calculable in 
terms of just one free parameter (the gravitino mass, m 3 / 2 ), 

2. the soft terms are real and flavor invariant, thus solving the SUSY flavor and 
CP problems, 

3. the soft terms are actually renormalization group invariant [IT], and can be 
calculated at any convenient scale choice, 
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4. the scale of the gravitino mass is too high to affect the Big Bang neucleosyn- 
thesis (BBN) bound and cosmological gravitino problem which was the main problem 
in SUGRA model. 

Supersymmetry breaking effects in the observable sector have in this framework 
a gravitational origin. In ordinary gravity-mediated supersymmetry breaking model, 
the supersymmetry breaking is transmitted from HS to OS via tree level exchanges 
with gravitational coupling. But in AMSB, the HS and the OS superfields are as- 
sumed to be located in two parallel but distinct 3-branes and the 3-branes are sep- 
arated by bulk distance which is of the order of compactification radius, r c . Thus 
any tree level exchange with mass higher than the inverse of r c is exponentially sup- 
pressed. So, the supersymmetry breaking is propagated from the HS to the OS via 
loop generated superconformal anomaly. The soft SUSY breaking terms related to 
gauginos and sleptons are calculated to be, 

Mi = ^1113 , (1) 
9i 2 

2 1 ( 9 1 a , d l \ 2 



A v = -^m 5 , (3) 

V y 2 

where M, is the gaugino mass term, tuq is slepton mass and A y is the trilinear 
parameter. 7 is the anomalous dimension and (3 is the beta function of this theory. 
7 and (3 are defined as, 

_ dlnZ _ dg = dy 
dlnpi ' 9 dlnji ' v dln/j, 

where Z is the renormalization constant for the gauge coupling, \i is the Higgsino 
mass, where (3 g and (3 y are, respectively, the gauge coupling and Yukawa coupling 
/3-functions, and their correspond anomalous dimensions are denoted by 7. Another 
feature of AMSB is that slepton mass-squared terms are negative giving to tachy- 
onic states as seen from equation. The problem is tackled by adding an universal 
mass-squared term Wq to all the squared scalar masses in the minimal extension to 
this theory, namely, minimal anomaly mediated supersymmetry breaking (mAMSB) 
model (O [13]). A sparticle spectrum in this model is fixed by three parameters, ms 
which is gravitino mass, tan/3 which is the ratio of the vacuum expectation values of 
the two Higgs fields (H® and H®) and sign(/i), where \x is the Higgsino mass. Also, 
an universal mass squared term (wiq) is needed to make all sparticle positive. So, 
ultimately, with mo, four parameters are needed to generate spectrum in mAMSB. 
So, we can generate various LSP neutralino masses out of these four parameters in 
this model. The neutralino is the lowest mass eigenstate of linear superposition of 
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the photino (7), zino (Z), and the two Higgsino states (Hf and H®) [14], written as, 

X = «i7 + CL2Z + a 3 H<? + a 4 H$ . (5) 

This state is obtained by diagonalising the neutralino mass matrix, which is given 
by, 

/ Mi —m z cos[5sin9 w m z sinf3sin8 w \ 

M 2 m z cosf3cos9 w —m z sinf3cos9 w 

—m z cosf3sin9 w m z cosf3cos9 w — ji 

\ m z sinj3sin6 w —m z sin(3cos9 w —ji o j 

in the basis (7 Z ffj 5 if 2 °) . 

Here /3 is the ratio of vacuum expectation values between the two Higgs doublets, 
m z is the mass of the Zq , 9 W is the weak mixing angle, and Mi ,/x are the U(l) 
and SU(2) gaugino and Higgsino mass parameters, respectively 

But there are certain phenomenological bounds on the parameters space, i.e, 

1. A lower limit on m? coming from the lower bound of m"? n = 86 GeV on 

2 x ± 

charginos causing direct search at the CERN-LEP. 

2. For a certain ms, there is a lower bound on mo below which f is LSP or 
sleptons are observables. 

3. For some choices of SUSY parameters, unbounded from below (UFB) directions 
of scalar potential are obtained and that parameter space region is not allowed. 

In a recent work Vasquez et al. [15] has given a detailed analysis of the allowed 
parameter space for a neutralino dark matter in the framework of NMSSM model. 
In their case the dark matter (neutralino) mass was within the range of ~ 80 GeV 
and hence the energies of the gamma rays from such dark matter annihilations can 
be probed by FermiLAT [16] experiment. In the present calculation, we instead 
consider the neutralino dark matter in mAMSB model mentioned above. Some of 
the earlier works on dark matter phenomenology in AMSB model include Baer et 
al. [lTj, Moroi et al. [IE], Ullio et al. PS] etc. In Refs. [T7] and [IS] the 7 flux 
from the galactic centre are discussed and although neutrinos from the neutralino 
annihilations are mentioned in Ref. [17] but they have not discussed elaborately. 
Moreover only two halo models are considered for their analysis. In an another 
earlier work ([20]), a neutralino dark matter in AMSB model is studied to obtain the 
region in scalar cross section (a scalar - parameter space. But in this case WMAP 
limit has not been taken into account. In Ref. [21], the 7 signal from galactic centre 
region due to dark matter annihilation is addressed mainly for the case of GLAST 
[22] satellite-borne experiment. Ref. [23] discusses the the 7- flux from galactic centre 
region, originated by dark matter annihilations. The authors made the analysis with 
different particle dark matter candidates with reference to MSSM, Kaluza-Klein extra 
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dimensional model etc. for different halo profiles and taking into account the Fermi- 
LAT experiment. But the neutrinos as dark matter annihilation products are not 
addressed. In another work by Allahverdi et al [23] considered MSSM and U(1)b-l 
extened MSSM model for dark matter candidate and calculated 7 and neutrino fluxes 
from galactic and extra-galactic origins by annihilating dark matter. But they have 
considered only NFW halo profile and they have not shown the neutrinos flux for 
different neutrino flavours. Moreover, no detailed comparison of their results with 
high energy neutrino or gamma search experiments. 

The phenomenology of the AMSB model (adopted in this work) was investigated 
by Datta et al. [25], among others, in an earlier work where they have calculated the 
allowed supersymmetric parameter space for different values of tan (3. In this work, we 
adopt their results for supersymmetric parameters such as m and ms and obtained 

neutralino masses upto TeV scale using the allowed region of m^-ms parameter space 

2 

shown by Datta et al. [25]. The relic densities for such dark matter were then 
computed using these SUSY parameters and they were compared with the WMAP 
results. The parameter, thus constrained further was then used to calculate the spin 
independent and spin dependent cross sections (o" sca tt) for different neutralino masses 
(m x ) (obtained using the restricted parameter space) in the cases of neutralino- 
nucleon scattering processes. These processes are essential for direct searches of dark 
matter. We found that the allowed m x — <T sca tt region, thus obtained, are found to 
be within the allowed limits of most of the direct detection experiment results. 

The indirect detection of dark matter involves the detecting the particles (and 
their subsequent decays) or photons produced due to dark matter annihilations. 
These annihilation products can be fermions or 7 photons. The dark matter parti- 
cles if trapped by the gravity of a massive body like sun or galactic centre they can 
annihilate there to produce these particles, in this work we investigated gamma rays 
from such annihilations of dark matter from the galactic centre and galactic halo re- 
gions. Using the constrained mAMSB parameter space discussed above we found the 
gamma ray flux from the galactic centre. We have also calculated the galactic neu- 
trino flux from such annihilations of mAMSB neutralino dark matter. These studies 
are performed for different galactic dark matter halo profiles. We found that the 
gamma spectrum from galactic centre and halo produced by neutralino dark matter 
within the framework of the present mAMSB model, is highly energetic. The exper- 
iment like HESS [2m [27], that can probe high energy gamma rays and which being 
in the southern hemisphere has better visibility of the galactic centre will be suitable 
to test the viability of the present dark matter candidate in mAMSB model. 

The possibility of detecting neutrinos from galactic halo from dark matter annihi- 
lations are also addressed with reference to Astronomy with a Neutrino Telescope and 
Abyss environmental RESearch (ANTARES) [25] under sea neutrino experiment. 
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The paer is organised as follows. In section 2 we discuss the calculation of relic 
densities of mAMSB neutralinos for the parameter space obtained from [25]. The 
relic densities are then compared with the WMAP results. The parameter space 
thus constrained further by WMAP is then used to calculate the spin dependent 
and spin independent scattering cross sections. They are compared with the existing 
direct detection experiment limits. These are discussed in section 3. In section 4 the 
indirect detection of the mAMSB dark matter from their annihilations at galactic 
centre and halo are discussed. To this end the gamma signals and neutrino signals 
are addressed. Finally in section 5 we give discussions and conclusions. 



2 Relic Abundance Calculation 

In order to calculate the relic abundance of the LSP, \i one needs to consider annihila- 
tion of N supersymmetric particles with masses rrtj (z=l,2,..,iV) and internal degrees 
of freedom gi respectively. If, R parity is conserved, then all the heavier supersym- 
metric paricles will decay to LSP. The relic abundance is obtained by numerically 
solving the Boltzmann's equation, 

diTi 

where n is the total number density of all the supersymmetric particles 

n = Ej7ij , 

and n eq is the value of n when the particles for dark matter candidate were in chemical 
equilibrium. At this epoch the temperature T of the universe was greater than Tf 
(T > Tf), the freeze out temperature of the particle considered. At a temperature 
below the freeze-out temperature Tf, the particles falls out of chemical and thermal 
equilibrium and their co-moving number density becomes fixed or "frozen". In Eq. 
El H denotes the Hubble parameter and (av) is the thermal average of the product of 
annihilation cross section and the relative velocity of the two annihilating particles. 

W 0') 

( av ) = ^((TijVtj) e ^ 2 eq , (7) 

with 



1,3 ec l 



(Pi-Pj) 2 - rnjm 2 - 
EiEj 



In the above, (pi,Pj) and (Ei, Ej) are the momenta and energies respectively for the 
ith and jth particles. The entropy conservation principle enables one to work with a 
useful quantity called abundance, Y = n/s [29] where s is the total entropy density 
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of the universe. Defining a dimensionless quantity x = m x /T, with m x being the 
mass of LSP, Eq. [6] can be written in the form 

In Eq. [HI Y eq is the value of Y when n = n eq . With Hubble parameter H = <J §7rGp, 
G being the gravitational constant, the total energy density (p) and the total entropy 
density (s) of the universe are given by [29] 

P = 9e ff (T)^T* (9) 

and a = /i e// (T)^T 3 . (10) 

In Eqs. M and [TU] 9eff, h e ff are the effective degrees of freedom for the energy and 
entropy densities respectively. Substituting Eqs. [91 Eqs. [10] and the expression for 
H in Eq. [HJ one obtains the evolution equation of Y as 



dx I 7T / X 2 N ' V 



g) ^-^M(^ 2 -0, (ii) 



where g l J 2 is defined as [29] 



91/2 = S + \£r& 1 • (12) 



The expression for Y eq is given by 

Y ^ - &£m*4 K > (?) ■ < 13 > 

where we sum over all supersymmetric particles denoted by % with mass rrii and 
internal degrees of freedom g^. K n (x) is the modified bessel function of the second 
kind of order n. Y eq (T) falls very rapidly as temperature decreases. The thermally 
averaged annihilation cross section, (crv), is calculated using a Maxwell-Boltzmann 
velocity distribution and given by, 

{aV) 2T{H i g i m 2 K 2 {m i /T)f ' ( ] 

where is the total cross section for annihilation of a pair of particles under con- 
sideration and pij representing the momentum of the incoming particles in their 
centre-of-mass frame and is given by, 

(s — (rrii + rrij) 2 )(s — (rrii — rrij) 2 ^ 



1 

Pij 



2 



(15) 
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and the summation runs over all particles. In Eqs. and [151 s denotes the relevant 
Mandelstam variable. The thermally-averaged cross section, (av) must include all 
channels by which \ can interact, including coannihilation with other particles, in 
which the number densities of both species are important. We also mention that the 
heavier particles contribute less in (av) due to strong Boltzmann suppression factor, 

fm i + m j -2m x \ 
B Ki exp I ^ -J (16) 

We have taken, B to be 10~ 5 to involve all coannihilation channels of heavier SUSY 
particles although in relic density calculations. 

Integrating Eq. ITT1 from x = xq = m/To to x = Xf = m/Tf, where To is the 
present photon temperature (2.726° K) we obtain Y (value of Y at T = T ) which 
is needed to compute the relic density. We have solved Eq. [TT] numerically. The 
following aprroximations used to solve Eq. [TTJ 

• 1. At small x (high T), the abundance of lightest SUSY particles (LSP) are 
almost in equilibrium and the temperature variation of the deviation from equi- 

y y 

librium abundance is negligible, i.e., Y ~ Y eq and — ^P- »s 0. Thus, the evolu- 
tion equation reduces to, 

dln(Y eq ) _ /45 \ ~ 1/2 gl /2 m x 



dx 



(-g\ 9 -^(av)Y eq 6(5 + 2) , (17) 

V 7T J x £ 



where 5 is some small constant coming from the definition of freeze-out tem- 
perature Tf. 

• 2. At temperature below Tf, equilibrium abundance, Y eq falls much below Y, 
as seen from Eq. [TTJ and can be neglected in the abundance evolution equation. 
Thus, Y is obtained from the relation, 

i i /45 \ _1/2 r° 9* /2 ( x 



r„ y, m * \~ G j l -ir-™** < 18) 

The relic density of LSP, in the units of critical density, p cr = 3H 2 /8nG, can be 
expressed as 

n x = w = Hb^o (19) 

Per Per 

where Sq is the present entropy density evaluated at T . Finally, knowing Y , we can 
compute the relic density of the dark matter candidate, from the relation [2Tj] , 

Q y h 2 = 2.755 x 10 8 -^-y . (20) 
x GeV 



S 



In the above h is the Hubble constant in 100 Km sec" ^Mpc -1 unit. The WMAP 
survey combining with recent observations of largescale structure provides the con- 
straints on the dark matter density Q DM h 2 as 

0.099 < n DM h 2 < 0.123 . (21) 

where Qdm is the ratio of dark matter density to the critical density p c = 1.88h 2 x 
10 _29 gcm~ 3 . 

In the present work we calculate the relic densities for the dark matter candidate 
neutralino in mAMSB model and compare our results with the WMAP bound. The 
allowed parameter space in the present SUSY model is thus extracted by WMAP 
results. 




0.0001 r 



1e-05 1 1 1 — 1 1 — 1 1 — 1 1 — 

1 10 100 1000 10000 

Mass of LSP neutralino, ( GeV) 

Figure 1: Scatter plot of Mass of the LSP neutralino, (m x ) vs. relic density (Qh 2 ). 
The cyan and pink line represents the WMAP upper and lower bounds on dark matter 
relic density respectively and the blue dotted zones corresponds to the mass range 
satisfying the WMAP limits. 

As mentioned earlier, the parameter space of the mAMSB model is defined by the 
four parameters, namely m 3 / 2 , m , tan/3 and sign(/x). The whole parameter space 
defined by the above parameters and constrained by the allowed region of mo — 7713/2 
as given in Ref [25] is used to calculate the relic density Q x (or Q x h 2 ) and the results 
are then compared with the WMAP results. 

The relic density in the present formalism of SUSY model is computed using the 
code micrOMEGAs [30]. We thus obtain the relic density for the scanned SUSY pa- 
rameter space discussed above. We find that the generated LSP neutralinos span 
very large range of mass. Each generated LSP neutralino mass gives rise to differ- 
ent annihilation cross section due to their annihilations to different standard model 
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particles and also co-annihilation processes. So, due to different (crv), the evolution 
of corresponding Boltzmann equation (Eq [8]) will be different giving rise to various 
relic abundances. 

In Fig. [U the variation of relic densities for different LSP neutralino masses are 
shown. The scatter plots in Fig. [T] correspond to the allowed parameter space of 
mo — ni3/2 obtained from (25]. The WMAP bound is superimposed on this scatter 
plot in Fig. Q] and the regions of agreement of the present calculational results with 
WMAP data are identified by blue coloured area in Fig. [TJ From Fig. (TJ we obtain 
two different neutralino mass regions satisfying the WMAP bound. One region is 
around 1 TeV and the other region at a somewhat higher range of ~ 2 TeV. 



WMAP relic lower bound 
WMAP relic upper bound 




Figure 2: Constraints on SUSY parameter space from WMAP limits in present SUSY 
model. The gaugino mass parameter niQ are shown by the colour index where mo 
varies from blue coloured region to yellow region as its mass increases. 

In order to elaborate how the WMAP bound constrains the SUSY parameter 
space in the present mAMSB model, we make a 3-D colour coded plot in Fig. [2J 
where the variation of relic density Qh 2 with the simultaneous variations of all three 
SUSY parameters namely tan /3 and mo are furnished. In Fig. [2J, the parameters 
m3/2 and tan/3 are plotted along X and Y axes respectively while the variation of 
gaugino mass m is shown in colour coded display whereby the colour reference deep 
blue denotes the lower value of m and increases towards the yellow zone in the plot. 
The corresponding variation of Qh 2 is shown alonf Z axis. The WMAP limits are 
shown in Fig. [2] by two meshes separated by the WMAP limit along Qh 2 axis. One 
observes from Fig. 2 that a very small region of the m 3 / 2 ~ m o ~ tan/3 parameter 
space is allowed by WMAP. Thus WMAP limit further constraints the m 3 / 2 — m 
parameter limits obtained from Ref. [23]. From Fig. [2Jit is also clear that only higher 
values of m (~ 10 - 12 TeV), and m 3/2 (~ 650 - 700) TeV could satisfy the WMAP 



10 



limits. 



Constrained by WMAP limits 



1e-25 
g 1e-26 



1e-28 1 ' 1 ' 1 1 

500 1000 1500 2000 2500 

( GeV) 

Figure 3: Plot showing the variation of annihilation cross section times relative ve- 
locity of annihilating neutralinos (o~v) with the mass of the LSP neutralino, (m x ). 
The green zones are the WMAP allowed regions 



In Fig. [31 we show how the annihilation cross sections vary with the neutaralino 
dark matter mass (m x ) in the present model. The WMAP allowed mass region is 
also shown by green colour. The av for the allowed zones (marked green) are seen 
to be around the value ~ 10~ 26 cm 3 sec _1 satisfying the "WIMP miracle" . 

The variations of freezeout temperatures (Tf) of LSP neutralino for the mass 
range obtained using the SUSY parameter space discussed earlier, are shown in the 
scatter plot of Fig. HI The neutralinos that satisfy WMAP relic density results are 
shown as blue in Fig. HI As in Fig. 1, in this case also one observes two such regions, 
one is around Tf ~ 80 — 86 GeV (more populated) and the other (fewer candidates) 
is at a lower region of around Tf ~ 40 GeV. 

3 Direct Detection 

The direct detection of dark is based on the principle that, the WIMP scatters off the 
target nucleus of the material of the detector causing the nucleus to recoil. The signal 
generated by the nuclear recoil (generally ~ keV) is measured for direct detection. In 
the direct detection experiments, attempts are made to give a bound in the m x — o~ sc 
space (m x being the mass of the dark matter and a sc is the dark matter-nucleus or 
dark matter- nucleon scattering cross sections). Different techniques are adopted by 
different direct detection experiments in order to measure the nuclear recoil energies. 
Some experiments that use Ge, Si or Nal as detector material use scintillation, phonon 
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Mass ofLSP neutralino, m y (GeV) 



Figure 4: Plot showing the variation of freezeout temperature (Tf) with the mass of 
the LSP neutralino, (m x ). The blue dotted zones are constrained by WMAP limits 
on dark matter relic 



or ionization techniques. In another class of detectors like Time Projection Chamber 
or TPC detectors, the drifting of ionized charges, produced by recoil nucleon of the 
detector gas (generally Xenon or argon), produce the track from which the direction 
of recoil can also be measured. Some of the ongoing direct detection experiments 
include DAMA [31] (Nal), CDMS [321 [33] ( 73 Ge), PICASSO [34] (CS 2 ), XENON 
[35j[36], COUPP [37] etc. They give different limits on scattering cross sections for 
different dark matter mass. 

For direct ditection case, we consider that the dark matter particles make an 
elastic scattering with the detector nucleus. The observation of such a process would 
require only a look for the nuclear recoil of detector material. 

The dark matter-nucleus interaction required to calculate the scattering cross 
sections, can be axial- vector (spin-dependent) or scalar (spin-independent). The 
target nucleus, with zero ground state spin gives rise to spin independent interaction. 
On the other hand, spin-dependent interactions are for the nuclei with unpaired 
nucleon that gives rise to non-zero ground state spin. Edelweiss [38], DAMA, CDMS 
SuperCDMS [39], XenonlO [35], XenonlOO [36], Zeplin [HI EH], Warp [42], KIMS [43], 
CoGeNT [44] are using detectors made of heavy nuclei (Ge or Xe) to search scalar 
interactions. On the other side, NAIAD [45J, SIMPLE [42], PICASSO, Tokyo/NaF 
[47] are using light nuclei to detect spin-dependent case. 

For non-relativistic limit and for Majorana fermionic WIMPs, the effective La- 
grangian for WIMP-nucleon elastic scattering is given by [48] 
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where A at, e and fcj's are different types of WIMP-nucleon couplings. At zero mo- 
mentum transfer limit, for the spin-independent case, the effective Lagrangian takes 
the form, 

L SI = Xn^^^niPn , (23) 

where suffix N stands for nucleon type. The squared amplitude averaged over the 
polarised states is given by, 

\A SI \ 2 = Q4(X N m x M N ) 2 , (24) 

Where m x , Mn are the dark matter particle mass and nucleon mass respectively. 
Summing over proton and neutron amplitudes of WIMP-neucleon scattering of a 
nucleus of mass number A and proton number Z, gives, 

\A SI \ 2 = U{m x M N f(\ p Z + X n (A - Z)) 2 . (25) 

Similarly, for spin-dependent case, the main contribution on Lagrangian comes from 
the last term of Eq. [22j The effective Lagrangian and the scattering amplitude for 
this case is given respectively by 

L S d = tNiixlvl^x^Nl^l^N , (26) 
\A SD \ 2 = 192(e N m x M N ) 2 . (27) 

For spin-dependent interaction at rest, the time component of the pseudovector cur- 
rent in the Lagrangian (Eq. 126]) can be neglected leaving a 3-dimensional spatial 
vector current proportional to the angular momentum J. The matrix element that 
is obtained summing over all spin states will also have a J dependence. The spin 
dependent WIMP-nucleus scattering amplitude A SD will also have a J dependence 
and is given by, 

|A SD | 2 = 256^±1(]T e t St) 2 (m x M N ) 2 . (28) 

i=p,n 

In the above, the factor Sf denotes the expectation value of the spin content of 
nucleon of type % and A is the mass number of the detecting material and this term is 
taken to be zero for even-even nuclei and Sf ~ 0.5 for odd nuclei. The spin-dependent 
and spin-independent cross sections for dark matter particle (x) with nucleon (N) 
are respectively given in compact forms as, 

Am 2 Mir 

a SD = X N x 3 | A SD|2 ^ ( 2 g) 



7r(m x + 



AmlM 2 N 
In) 



= -"V"^ x | A SI|2 (3Q) 
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Finally the spin independent and spin dependent cross sections for WIMP scattering 
to a pointlike nucleus at rest [7] are given by the following compact forms 



16m 2 x M 2 ^J A + 1 ( y- r , 

i=p,n 



'^^^T^^ 1 ' (31) 



Am 2 M 2 

^ = , * " x (X P Z + \ n {A - Z)f . (32) 
ii(m y + Mn) 



-C2 



spin-independent scattering cross-section satisfying WMAP bound 




Mass ofLSP neutralino, (GeV) 

Figure 5: The plot for Mass of the LSP neutralino, (m x ) vs. Spin Independent 
scattering cross section (cr ). 



We have computed both spin-independent and spin-dependent scattering cross 
sections of neutralino dark matter for a wide range of mass in this model respecting 
the mo — rriz/2 bound [25] . As the nucleon consists of both protons and neutrons, 
the WIMPs can be scattered off both nucleons. The contribution of loop diagrams 
alongwith the tree level diagrams has also been included for calculations of scatter- 
ing amplitudes for both SI and SD cases of \ — N scattering. These scattering cross 
sections for different neutralino masses are computed using micrOMEGAs [30] com- 
puter code. The results for both SI and SD cases are shown in Fig.()5]) and Fig(|6]) 
respectively as scattered plots for scattering with protons. 

The green scattered plots in Fig. give the spin independent scattering cross 
section a SI for various neutralino masses m x generated in the present AMSB model 
with the bound on mQ— m 3 / 2 parameter space given by Ref. [25]. The blue scattered 
plots in Fig. [61 on the other hand, are for the spin dependent case. The mass region(s) 
in this model that satisfy the WMAP results for relic density are superimposed over 
these two figures in order to constrain the m x — <t SI//sd space obtained from Figs. El 
|6j The blue patches in Fig. [5] and the pink patches in Fig. [6] represent the mass 
regions that satisfy WMAP results. Clearly, there are two different zones allowed 
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Figure 6: The plot for Mass of the LSP neutralino, (m x ) vs. Spin Dependent scat- 
tering cross section (a SD '). 



by the WMAP limits as expected from the discussions in Sect. 2. For the WMAP 
allowed lower mass region (around 1 TeV), the SI cross section (Fig. \5§ extends 
between ~ 1(T 9 - ~ 1(T 7 pb. The WMAP allowed higher mass region (around 2 
TeV) which spans larger region in m x — o~ sl space than the WMAP allowed lower 
mass region, is confined within SI cross section limit ~ 10~ n — ~ 10~ 14 pb in Fig. 
The pink regions in Fig(jS]) signify the WMAP allowed region. Here, the value of SD 
scattering cross section is coming to be higher than that of SI as it is expected from 
the theoretical perspective. 
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DATA listed top to bottom on plot 

CDMS I (SUF), 2010, shallow site results of Ge k Si detectors, SI 

— CoGeNT,2010,SI 

— CDMS E (Soudan), 2011, reanalyzed data from Oct06-Sep08, Ge detector, 2keV 

- - - NAIAD, 2005, final result, SI 
KIMS. 2007, 3409 kg-days Csl, SI 

CRESST 1, 2004, lOJkg-days. CaW04 detector, SI 
| DAM, 2000, NaI-0 to4 combined, 58,0O0ks-days. 3sisma, SI 

- - - WARP, 2007, 21, %i ke-days 40 keV threshold. 5 events 40-54 keV assume 

ZEPLFN ID, 2009, first science run, 847kg-days, SI 

Edelweiss II, 2011, Final Results, Ms-days, SI 
I I ZEPLIN III, 201 1, second science run, l344kg-davs. SI 

CDMS E (Soudan), 2004 to 2009 combined all data from Soudan, SI 

I I XENON10. 2007. Net 136ka-d 

XENON100, 2011, 100.9 live days of data, SI 

I I this work 



Figure 7: Limits on spin independent scattering cross sections set by various experi- 
ments 



Similarly, in Fig. |6]the WMAP allowed lower mass region (around 1 TeV) con- 
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strain the spin dependent cross section a limits in the range ~ 10 — ~ 10 
pb and for the region of around 2 TeV <r SD lies between ~ 10~ 10 to ~ 10 -9 pb. 
We mentioned in passing that we obtained similar nature for WIMP '-neutron elastic 
scattering. 




10" 



10 

WIMP Mass [GeV/c 2 ] 



DATA listed ton to bottom on plot 

ZEPLIN E 2009, first science ran. Bonn-CD nuclear potential. SD-proton 

■ ■ ■ Edelweiss. 2005. 4.8 kg-days. 73-Ge detector. SD-proton 
ZEPLIN 1 2007. fetalis, 31.2 live days. SD-proton 

x x x CRESST 1. 2004. (est). 10l£-days. CaWo4 detector. SD-proton 

■ ■ ■ CDMS II (Soudan). 201 1. reanalyzed data from OctOfrSepOl Ge detector. 2keV 



SIMPLE. 2005. usine four superheated droplet detectors with a total exposure of 
| DAMA. 2003 est. Nal detector. SD-proton 
■ ■ ■ XENON10. 2008. 58.6 Live days. SD-proton 
NAIAD. 2005. fmal result. SD-proton 
xxx KM. 2007. 3405 kg-days Csl. SD-proton 
PICASSO, 2010, 13.75±0.48kg— days. SD-proton 
COUPP. 201L using 3.5kg CF^3I detector. SD-proton 
xxx^ SIMPLE, 201 1, Final analysis of Phase n, SD-proton 



Figure 8: Limits on spin dependent scattering cross sections set by various experi- 
ments 



Figs. [7]and|8]show respectively various upper limits in dark matter mass - SI or SD 
scattering cross section (m x — cr SIorSD ) plane set by different ongoing direct detection 
experiments. The WMAP-allowed regions from the present model for neutralino 
dark matter are superimposed on them for comparison. The experimental limits 
are obtained from the compilation given in Ref. [19]. The names of the different 
experiments are furnished as legends in the figures. It is obvious from Figs. [7] and 
E] that the allowed parameter space for the considered AMSB model is within the 
allowed limits of the experimental bounds. 

4 Indirect Detection 

Weakly interacting dark matter in our galaxy can be trapped inside massive heavenly 
bodies like galactic centre or the sun due to the gravity of these bodies. The dark 
matter particles, in course of their passage through such massive bodies undergo 
elastic scattering off the nuclei present there as a result of which their velocity deplete. 
If their velocities fall below their escape velocities from such massive objects, the dark 
matter particles are trapped. These trapped dark matter particle (x) may undergo 
the process of pair-annihilation producing primarily b, c and t quarks, r leptons, 
gauge bosons, etc. (xx ~~ >" qq,l + l~ .uu, ZZ,W + W~ , ...). The annihilation products 
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depend on the mass and composition of the dark matter. Neutrinos and antineutrinos 
can be produced by the decay of primary annihilation products or through direct 
annihilation. 

The main principle of indirect detection of dark matter is to detect and measure 
the fluxes of standard model particles produced from the annihilation of dark matter 
trapped by the gravitation of massive heavenly bodies. 

Recently many new results from indirect DM searches have been released. An 
interpretation of these excesses related to astrophysical processes from any galactic 
or extragalactic sources is still not very clear. The products from the annihilation of 
dark matter particles in massive bodies such as in galactic centre may explain such 
excess signals. 

There are a lot of satellite borne experiments that look for gamma rays or anti- 
matters in cosmos. Some terrestrial experiments are also suited for looking at cosmic 
gamma rays, neutrinos etc. 

Such experiments include PAMELA [50] that confirms an excess in positron frac- 
tion in agreement with earlier indications by HEAT [51] and AMS01 [52]. Other 
satellite borne experiments like FERMI [53] and ATIC [51] report an excess in total 
electron and positron spectrum at energies of several hundreds of GeV's, much higher 
than that of PAMELA search. The cosmic gamma rays from the galactic sources and 
from galactic centre are measured in a wide range of energies by INTEGRAL (< ~1 
MeV) [55], EGRET [5§], FERMI [57], GLAST (-10 MeV- -100 GeV) [22], HESS 
HUES], MAGIC [5B], Whipple/Veritas [53], CANGAROO (> -100 GeV) [§0] etc. 

In this work we mainly focus on the gamma ray and neutrinos from dark matter 
annihilations at and around galactic centre (GC) as these regions have higher dark 
matter density and hence a promising site for the study of indirect detection of dark 
matter. 

Although GC seems to be the most obvious target, it is also one of the most diffi- 
cult areas to work with because of the complex and poorly-understood backgrounds 
[6Tj 152"] . for signals from around GC. and uncertain dark matter profile [63] IM] . 

The galactic gravitational potential leads to a higher dark matter density at the 
centre of Milky Way. The expected flux from the galactic centre depends on the 
distribution of dark matter in the galaxy. The dark matter density profile p(r) is 
assumed to be spherically symmetric. The differential flux of the outgoing particle 
of type i is given by 

P(E, 8) = ^ = Y -^ d -^(E)m Ail) (33) 
v ' ' dE ^ 8Traml dE y ' v ' ' y ' 

where the factor, a is 1 or 2 depending on whether the assumed WIMPs are self- 
conjugated or not respectively. It is also to be mentioned that the effect of this factor, 
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a on the above differential flux is much less significant in comparison to the dark 
matter density fluctuations in the innermost regions of Milky Way. Here we consider 
a to be unity as the neutralinos from the mAMSB model (the dark matter candidate 
chosen in the present work) are self-conjugated. In Eq. [331 a is the annihilation 
cross section of dark matter and v denotes the relative velocity of the dark matter 
particles. The quantity *%r(E) in Eq. [33] is the energy spectrum of particle % and 
J(0, Aft) is given by, 



j (e, An) 



An 



(p 2 (r(f,9)))df 



(34) 



line of sight 



With 9 being the angle subtended by the line of sight of an observer on the earth 
(along the length f) on R & — the distance between GC and the terrestrial observer 
(in solar system). The source to observer distance f can be calulated as 



2rR Q cos 



(35) 



In the above the dark matter annihilation region is considered to be at a distance 
r from galactic centre. Here we also mention that the galactic centre is assumed to 
be coincident with the halo centre). The solar system's position in the halo from 
the GC is given by Rq = 8.0 kpc. In Eq. [3H AQ is the solid angle over which the 
observation is to be made and p(r) is the dark matter density at a distance r from 
GC. Clearly the integration on the RHS of Eq. [34] is along the line of sight. Thus 
the astrophysical factor J in Eq. [321 has only a 9 dependence (alongwith AQ) and 
thus the differential flux J 7 can be expressed in terms of the angle 9 corresponding 
to different positions of the source in galactic halo with respect to GC. 

The dark matter density p(r) is related to the spherically symmetric halo profile 
of galactic dark matter by the equation 



p(r) 



Po^haloM 



(36) 



where po is the dark matter density at the galactic centre assumed to be 0.3 GeV/cm 3 
and -Fhaio( r ) is the halo profile of the galactic dark matter which can be expressed in 
a parametric form, 



F^\oir) 



Rc. 



■0 



Re 

a 



(37) 



In the above, a, a scale parameter and the other parameters a, /3, 7 take different 
values for different halo models which follow the above parametric form for -Fh a i . For 
example, for NFW halo profile [OS], a = 1, /3 = 3, 7 = 1 and a = 20 kpc, whereas 
the parameter set a = 2, (3 = 2, 7 = and a = 4 kpc represents isothermal profile 
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with core [67]. Again for the Moore profile [68] we have, a = 1.5, f3 = 3, 7 = 1.5 and 
a = 28 kpc. In Einasto halo profile [69] however, a different kind of parametric form 
is adopted which is given by, 



F^(r) = exp 



a 




where a is the parameter. In this work a = 0.17 is adopted. In what follows the four 
profiles are referred to as NFW, Isothermal, Moore and Einasto respectively. The 
galactic halo densities for these four halo models are shown in Fig. In the present 
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Figure 9: The variation of galactic halo density with radial distance for various halo 
models 



work we show the gamma ray and neutrino flux from galactic centre region for each 
of the four halo profiles mentioned above. 

4.1 Gamma Ray Flux Results 

There are generally two kinds of 7-ray emission from DM annihilation. In the first 
category 7 is produced directly from the annihilation final state particles which is 
called primary emission in which final charged leptons emit gamma ray or tt° decays 
to gamma ray after hadronization. The other kind is called secondary emission in 
which gamma rays are produced by interactions of final state particles with external 
medium or radiation field such as the inverse Compton effects etc. Here we consider 
only the first type of emission for which the relation (I3"3"|) holds. Annihilations can 
proceed directly into two photons [70j [71] or into a Z boson and a photon [72] via 
loop induced processes and the annihilation spectrum for both of these processes are 
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respectively given as, 



dN2 n 



dE 



dE 



(E) = 25{E 
(E) = s(e 



-m Y + 



m x ) 




) 



(40) 



(39) 



These processes give rise to monochromatic gamma ray lines. Here we calculate the 
gamma ray flux from the galactic centre as also from other places in galactic dark 
matter halo away from the GC. As discussed in the previous section, the positions 
away from the galactic centre are characterised by changing only the angle 9. This 
angle 9 in fact denotes the angle of sight from the observer with respect to the line 
of sight when the observer is looking directly at the galactic centre. The polarisation 
effect of final state gauge bosons (W^ and Z) and also the photon radiation effect 
which strongly affect the gamma ray spectra are also taken into account in the present 
work. The calculations are made for each the four halo profiles, referred to as NFW, 
Isothermal, Moore and Einasto and the results are furnished in the four figures, Figs. 
ITOk - [TOt i respectively. 

In Figs. fTOh - HUH , we plot the quantity E 2 x ^ for different values E, the 
energy of the emitted 7 rays from dark matter annihilations. We show the results 
for the cases when 9 = (galactic centre), 9 = 30° and 9 = 60° and are shown as 
red, green and blue regions respectively. One notices in Figs. [TOk - [TUb l that 7 flux 
for any particular value of the angle 9 are given as a pair of plots designated by 
the specific colour code (red, blue or green), assumed for the results corresponding 
to that particular 9. This is due to the fact that the WMAP data constrain the 
supersymmetric parameter space considered in this work, in two distinct zones as 
shown in Figs. [JJ [2] and each of the plots in every such pair of 7 flux in Figs. [TUk - 
[TOl d correspond to each of the WMAP allowed regions for the present AMSB model 
for cold dark matter candidates. It is clear from Figs. [TUh - [TUb l that calculations 
with different halo profiles yield different results for 7 flux. It is also to be noted that 
the flux from the galactic centre (9 = 0) is larger than the flux from other positions 
in galactic halo (corresponding to different values of 9 7^ 0) for each of the four halo 
models considered. 

The 7 fluxes are found to be almost of the similar order for the cases when 9 = 30° 
and 9 = 60° in each of the Figs. [TUh - [TUH . This reflects the fact that the DM halo 
profiles are almost flat in those regions. The Einasto profile has a finite (zero) central 
slope, unlike the NFW profile which has a divergent (infinite) central density. As it 
is not yet known which model provides the best description of the central densities 
of simulated dark- matter halos, we have taken these known models into account. 

The obtained 7-flux due to dark matter annihilation at the galactic centre for 
different halo model is compared with the observational results of The High Energy 
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Figure 10: Plot showing the variation of gamma ray flux with energies from the 
annihilation of dark matter for different galactic DM halo models and for different 
angles of sight, 9. The red lines describe the flux observed at 9 = 0°, i.e., from the 
galactic centre and the green and blue coloured regions are for the observations at 
9 = 30° and 60° respectively. The subfigures are for different commonly used dark 
matter halo profiles implemented in this work, a) NFW profile b) Isothermal profile 
with core c) Moore profile d) Einasto profile 



Stereoscopic System (HESS) experiment. Located in Namibia, the HESS experiment 
is designed to investigate high energy cosmic gamma rays (~ 100 GeV - TeV energy 
ranges) and it can also investigate the 7-rays in its observable energy range which 
can be due to the annihilation of cold dark matter particles. The results are given 
in Fig. [11] and Fig. [12j In each of the figures, the calculated flux are shown by two 
diferent regions corresponding to WMAP constrained two zones of dark matter mass 
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in the present mAMSB model (discussed earlier). It has been argued by Prada et 
al. in Ref. [73] that due to the infall of baryons at the galactic centre, the expected 
7 signal from dark matter annihilation at galactic centre will be boosted in case the 
dark matter consists of supersymmetric particles. In fact, considering neutralino in 
minimal supergravity (mSUGRA) model as the candidate for dark matter and with 
the NFW dark matter halo profile, they have demonstrated that the said boost can 
be of the order of 1000. In Fig. ITTk . the 7 flux from the galactic centre as calculated 
from the annihilation of neutralino dark matter in present mAMSB model assuming 
the NFW profile, is compared with the HESS results. The solid angle at which the 
HESS experiment looks at the galactic centre is ~ 10~ 5 sr, a value which is also 
adopted in the present calculations to obtain the results shown in Fig. [H] and Fig. 
[T2l The two regions in It is evident from Fig. ITTa that the 7-flux obtained from the 
calculations is much less than the HESS results for the energy range given by the 
model with WMAP constraints. In Fig. [TTb we show a representative plot where the 
calculated 7-flux is multiplied ("boosted") by a factor of 1000 and then compared 
with the HESS results. Fig. [TTb shows that the "boosted" flux is in the similar ball 
park of HESS results which seems to satisfy the claim made in Ref. [73] • In Figs. 
[T2"a . IT2"b and [T2"c . we show similar comparisons with HESS results for calculations 
made with Einasto, Moore and isothermal halo profiles respectively. One observes 
from Fig. [12] that both for isothermal and Einasto profiles, the calculated fluxes 
are below the HESS results while for Moore profile, they are comparable with HESS 
results. Both NFW and Moore profiles are cuspy in nature and they essentially differ 
by the values of the parameters a, (3, 7. On the other hand both the Einasto and 
isothermal profiles are non-cuspy in nature while the latter is a flat halo profile. One 
needs to increase the calculated flux by a factor ~ 10 2 for the former case while the 
calculated flux for the isothermal profile needs a boost of ~ 10 5 to be in the regime 
of HESS observational results. 

4.2 Neutrino Flux Results 

As discussed earlier in this section neutrinos can also be produced by the annihilation 
of two neutralinos - the present dark matter candidate. These trapped dark matter 
at the galactic centre produces primarily b, c, t quarks, r leptons, gauge bosons, etc. 
through the process of pair-annihilations. The neutrinos can be obtained from the 
decay or pair annihilation of the primary products. The neutrinos can also be pro- 
duced directly from the annihilation of two AMSB neutralinos (xx y v) mediated 
by Z, sneutrino [y] etc. In this work we investigate the muon neutrino (z/ M ) flux from 
the galactic centre due to the annihilation of such neutralinos in the prsent AMSB 
model and its possible detection prospect in an earthbound detector. Searches for 
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Figure 11: Plot of gamma flux with the energy for dark matter annihilation at the 
galactic centre and comparison with the HESS experimental data for NFW profile a) 
without baryonic compression and b) with the baryonic compression and ~ 10 3 flux 
enhancement 



neutralino annihilation into neutrinos is subject to extensive experimental investiga- 
tions in view of the neutrino telescopes like IceCube [71], Baikal [73], NESTOR [75] . 
ANTARES [28J. The calculation of flux of neutrinos coming from GC are similar to 
that of gamma rays as both are electromagnetically neutral particles. So, they are 
not affected by the irregularities of galactic magnetic fields or any magnetic turbu- 
lances. Also, they do not suffer any energy losses from inverse compton effect or from 
synchrotron radiation. For the present case we calcuate the possible muon (/x) signal 
from these neutrinos at ANTARES neutrino telescope [28] installed in the sea-bed 
off France coast. 

Firstly, we calculate the neutrino flux at the galactic centre for all the four halo 
models considered. The neutrino flux for the halo models can be obtained using 
similar equations (Eqs. I3"3l - I3"7j) that is used for obtaining 7- flux. The z/-flux for each 
of the three flavours namely u e , and v T are calculated separately for three values 
of the angle 9 (Eq. [33] and discussions earlier in this section) namely 9 = 0°, 9 = 
30°, 9 = 60°. The results are furnished in Figs. [T3] - E3 The three figures Fig. [121 
Fig. HU and Fig. [TS] correspond to 9 = 0°, 30° and 60° respectively. In these figures 
figures the v e flux , the flux and the v T flux are shown by red, blue, green colour 
in Fig. [131 pink, yellow and turquiose colour labels in Fig. [14] and black, yellow and 
orange in Fig. [15] respectively. The two flux regions for each of the neutrino flavours 
are for the two different allowed dark matter mass zones in this model obtained from 
WMAP results. The v flux for different dark matter profiles considered here, exhibit 
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Figure 12: The plot of energy vs. gamma flux from dark matter annihilation at the 
galactic centre and compared them with the HESS experimental data for a) Moore 
profile, b) Isothermal profile with core and c) Einasto profile 

similar trends as for the case of 7 flux in the sense that the flux is more for Moore 
profile and gradually decreases for Einasto, NFW and isothermal profiles. 

The neutrinos, while reaching the earth from the galactic centre will undergo 
flavour oscillations, whereby the flux of a particular flavour, say u^, will be modified 
on reaching the earth from the galactic centre. Since the baseline length L is very 
large in this case in comparison to oscillation length, the osciillation part is averaged 
out. Thus in the limit L — > 00, the probability that a neutrino with flavour a will 
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Figure 13: neutrino flux for three flavours (i> e , and v T ) for different energies from 
the annihilation of dark matter at from the galactic centre. The red, blue and green 
patches describe the fluxes cooresponding to v e , and v T respectively 



oscillate to flavour (3 is given by 

P(l/ a -> Vp] L = OO) = S a(3 - ^ U ai U fSiUajU^ 

= \U m \ 2 \U Pl \ 2 , (41) 

where a, (3 denote the e, \i or r, the different flavour indices and i,j = 1, 2, 3 denote 
the mass indices of three neutrinos. In the above, the oscillation part (~ Am 2 j(L/ E)) 
is averaged over due to large L/E. The mass- flavour mixing matrix U is denoted by 
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Figure 14: neutrino flux of three flavours (v e , and v T ) for various energies for 
angle of sight, 9 = 30° from the galactic centre. The pink, yellow and cyan coloured 
zones describe the fluxes cooresponding to v e , and v T respectively 



and 



U 



( u el u e2 u e3 \ 

Ufii Ufj2 U^s 



S13 



(43) 



V U Tl U r2 U r3 J 
In fact U is the usual MNS mixing matrix given by 

/ C12C13 S12S13 

U = ~~ S12C23 ~ C12S23S13 C12C23 — S12S23S13 S23C13 . (44) 

\ s 12 s 23 — C12C23S13 — C12S23 — Si2C23S 13 C23C13 J 

In the above, s, c denote sin 6, cos 9 respectively and 6 I 12 , #23 and #13 are three mixing 
angles for three neutrino species. We consider here no CP violation in neutrino sector. 
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Figure 15: neutrino flux of three flavours (v e , and v T ) for different energies for 
angle of sight, 9 = 60° from the galactic centre. The black, yellow and orange regions 
describe the fluxes cooresponding to v e , and v T respectively 



From Eq. [41] probability P can be written as 

P = XX T 

where the matrix X is given by 
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Hence, the oscillated flux of the neutrinos (of three flavours) at the detector is given 
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where the quantities in the RHS with superfix denote the initial neutrino fluxes. 

In the present work we estimate the muon yield for such a flux from galactic 
centre at ANTARES neutrino detector. ANTARES is a deep sea neutrino telescope 
and is basically a water Cerenkov detector, which detect the neutrinos by detecting 
the Cerenkov light of a charged lepton that is produced by the charged current 
scattering of neutrino off the sea water. The telescope consists of several vertical 
strings of around 350 metres long, each of which is fixed with 75 optical modules 
containing photomultiplier tubes. The strings are installed at the Mediterranian 
sea bed at a depth of around 2.5 Km off the French coast of Toulon. Designed to 
detect neutrinos with high energy (~ 100 GeV to ~ 100 TeV) of generally cosmic 
origin, this telescope looks in the direction of southern hemisphere. In fact due 
its position, ANTARES is very much suitable for observing the galactic plane and 
galactic centre. In the present context, the dark matter mass from AMSB model 
that is allowed by WMAP, is in the region of ~ 1 TeV - ~ 2 TeV and since neutrinos 
from the annihilation of such dark matter at galactic centre is being studied, this 
telescope is best suited for the purpose. 

The spectrum of muon yield, Q^(E Uix ), for different enrgies E v at ANTARES can 
be estimated using the relation 

= 0vW(^,)> ( 48 ) 

where, A e g jU is the neutrino effective area for ANTARES telescope and is obtained 
from Ref. [28]. 

The U/j, flux, (j) v , at the earth from the galactic centre is calculated using Eqs. 
HH - H7] with <P° v ^i the flux at the source, given in Fig. [T3] for different halo profiles. 
Note that, at galactic centre is considered only for the case 6 = (see earlier in 
this section) in the present calculation. The values of three neutrino mixing angles 
in Eq. IHare taken to be 9 12 = 34.0° [77], 9 23 = 46.1° [77J and 9 13 = 9.2° [78]. The 
results for estimated yield of muon spectrum $^(.2^ ) at ANTARES is shown in Fig. 
[TBI for all the four halo profiles considered. The estimates are shown for 5 year run 
of the telescope. It is seen from Fig. [TSJ that while NFW profile predicts very large 
yield the same using isothermal profile is rather low. The NFW profile has a cuspy 
structure whereas the isothermal profile gives a flat halo. 
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Figure 16: Estimated /i events for five year run at ANTARES neutrino telescope for 
different v„ energies obtained from dark matter annihilations at the galactic centre. 



5 Summary and Conclusion 

In this work, we have investigated the phenomenological implications of dark matter 
coming from a very well known SUSY breaking model, namely minimal anomaly 
mediated supersymmetry breaking (mAMSB) model. The suitable candidate in this 
model is the neutralino stabilized by the conservation of R-parity in SUSY theory. 
We have randomly scanned the parameter space of this model within the theoretical 
bounds (given in ref. [25]) of the parameter space of this model and for each point in 
paremeter space, we obtain a neutral stable candidate (neutralino) of dark matter. 
We have found the relic densities of each of these candidates and hence obtained 
the mass zone(s) which satisfy the WMAP constraints. The allowed parameters 
determined from such constraints are then used to study the direct and indirect 
detections of the proposed dark matter candidate, neutralino, in the present mAMSB 
model. 

The scattering cross sections for the dark matter particles scattered off the nucleus 
of the detecting material are determined by the nuclear form factors and the dark 
matter-nucleon coupling. The two types of cross sections, namely spin independent 
(zero nuclear spin at the ground state) and spin dependent, are determined by the 
different form factors and dark matter-nucleon couplings. We calculate both the 
spin dependent and spin independent scattering cross sections with the constrained 
zone(s) of the present neutralino dark matter parameter space and hence compared 
our results with several recent ongoing direct detection experimental results. The 
calculated mass-cross sections, thus calculated, are found to be below the upper 
limits of many well known experiments. From experimental point of view, the future 
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advanced direct detection techniques may probe those regions in mass-cross section 
plane given from the model. We also studied studies in detail the possibilities of 
indirect detection of these dark matter candidates. 

We have computed the the gamma ray flux coming from the annihilation of dark 
matter in mAMSB model in the galactic halo. For this reason, we have taken several 
well known theoretically motivated dark matter halo profiles such as NFW profile, 
Moore profile, isothermal profiles with core and Einasto profiles and calculate the 
flux coming from different positions of the halo plane. As the allowed mass of the 
neutralino (dark matter) is high (~ GeV to ~ 1000 TeV), so the energies of the gamma 
rays from dark matter annihilations are also of that order. Therefore high energy 
gamma ray search experiments may verify the present model. For this purpose, we 
have chosen the HESS experiment and compared our results for different halo profile 
considered, with the observed 7-flux of this experiment. We find that the 7-flux for 
non-cuspy profiles like isothermal profile (flat) and Einasto profile, is orders below 
the HESS results whereas the cuspy profiles like Moore profile overestimate the HESS 
result. Calculations using the other cuspy profile, namely the NFW profile requires 
a boost of ~ 10 3 for comparison with HESS results. The Moore profile has a ^5 
dependence while NFW profile varies as -. Thus the former is steeper than the latter. 
The two different cuspy nature appears to influence the result. The other two profiles 
namely the isothermal (flat) profile and the Einasto profile yield the calculated flux 
for the present dark matter candidate to be below the HESS results. It is still a 
matter of investigation to understand whether the galactic centre has a flat profile or 
a steeper halo profile. The present analysis, within the framework of mAMSB model 
considered in this work for dark matter candidate, seems to suggest that the cuspy 
nature of the profile appears to explain the HESS data better than the flat ones. We 
also like to add that we performed similar calculations with another flat halo profile 
namely Burkert profile [801 E], but the calulated 7 flux is found to be even below 
than what is obtained for isothermal profile. 

Different flavours of neutrinos from the dark matter annihilation at galactic centre 
are also addressed in the present work. The flux and detection of muon species of such 
neutrinos are calculated for the neutralino dark matter in mAMSB model. Given the 
masses of such dark matter candidates the energies of such neutrinos will also be in 
the range GeV to TeV. The location of the galactic center with respect to earth is 
downwards. The high energetic muon neutrinos may produce muon by the charged 
current scattering off ice or water and may be detected by their Cerenkov lights. We 
have calculated the fluxes of neutrinos of different flavours for annihilations of dark 
matter in present model at galactic centre as also at the other two chosen positions 
in the halo. The results are shown for the four halo profiles considered. In order to 
estimate the detection yield of such neutrinos in a terrestrial neutrino observatory, 
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we have chosen the ANTARES under sea detector and calculated the muon yield 
for muon neutrinos from galactic centre in this model for all the four halo models 
considered. The calculations of neutrinos in case of different halo profiles also exhibit 
similar trend as those for the calculation of 7 flux. 

The value of thermal average of the squared halo density, (p 2 {r)) is generally 
greater than ((p(V))) 2 due to the influence of a probable clumpy structure of dark 
matter halo profile, F c (v), which is related to dark matter halo profile by, 

(P 2 (r)} = p 2 FU(r)F c (r) (49) 

The clump structure of dark matter halo gives enhancement factor. In the present 
study of different models of galactic halo structures, we did not consider any clumpy 
halo of dark matter. This study is for posterity. 

The WMAP allowed zone(s) for the mAMSB model for dark matter, are around 
(~ 1 TeV and ~ 2 TeV) which are high in mass regime like Kaluza-Klein dark matter. 
The future collider experiment may verify their existence. 
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